Detection and Molecular Characterization of the SARS-CoV-2 Delta Variant and the Specific Immune Response in Companion Animals in Switzerland.

In human beings, there are five reported variants of concern of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2). However, in contrast to human beings, descriptions of infections of animals with specific variants are still rare. The aim of this study is to systematically investigate SARS-CoV-2 infections in companion animals in close contact with SARS-CoV-2-positive owners ("COVID-19 households") with a focus on the Delta variant. Samples, obtained from companion animals and their owners were analyzed using a real-time reverse transcriptase-polymerase chain reaction (RT-qPCR) and next-generation sequencing (NGS). Animals were also tested for antibodies and neutralizing activity against SARS-CoV-2. Eleven cats and three dogs in nine COVID-19-positive households were RT-qPCR and/or serologically positive for the SARS-CoV-2 Delta variant. For seven animals, the genetic sequence could be determined. The animals were infected by one of the pangolin lineages B.1.617.2, AY.4, AY.43 and AY.129 and between zero and three single-nucleotide polymorphisms (SNPs) were detected between the viral genomes of animals and their owners, indicating within-household transmission between animal and owner and in multi-pet households also between the animals. NGS data identified SNPs that occur at a higher frequency in the viral sequences of companion animals than in viral sequences of humans, as well as SNPs, which were exclusively found in the animals investigated in the current study and not in their owners. In conclusion, our study is the first to describe the SARS-CoV-2 Delta variant transmission to animals in Switzerland and provides the first-ever description of Delta-variant pangolin lineages AY.129 and AY.4 in animals. Our results reinforce the need of a One Health approach in the monitoring of SARS-CoV-2 in animals.

Large-scale analysis of SARS-CoV-2 synonymous mutations reveals the adaptation to the human codon usage during the virus evolution.

Many large national and transnational studies have been dedicated to the analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, most of which focused on missense and nonsense mutations. However, approximately 30 per cent of the SARS-CoV-2 variants are synonymous, therefore changing the target codon without affecting the corresponding protein sequence. By performing a large-scale analysis of sequencing data generated from almost 400,000 SARS-CoV-2 samples, we show that silent mutations increasing the similarity of viral codons to the human ones tend to fixate in the viral genome overtime. This indicates that SARS-CoV-2 codon usage is adapting to the human host, likely improving its effectiveness in using the human aminoacyl-tRNA set through the accumulation of deceitfully neutral silent mutations. One-Sentence Summary. Synonymous SARS-CoV-2 mutations related to the activity of different mutational processes may positively impact viral evolution by increasing its adaptation to the human codon usage.

Unsupervised explainable AI for molecular evolutionary study of forty thousand SARS-CoV-2 genomes.

BACKGROUND: Unsupervised AI (artificial intelligence) can obtain novel knowledge from big data without particular models or prior knowledge and is highly desirable for unveiling hidden features in big data. SARS-CoV-2 poses a serious threat to public health and one important issue in characterizing this fast-evolving virus is to elucidate various aspects of their genome sequence changes. We previously established unsupervised AI, a BLSOM (batch-learning SOM), which can analyze five million genomic sequences simultaneously. The present study applied the BLSOM to the oligonucleotide compositions of forty thousand SARS-CoV-2 genomes. RESULTS: While only the oligonucleotide composition was given, the obtained clusters of genomes corresponded primarily to known main clades and internal divisions in the main clades. Since the BLSOM is explainable AI, it reveals which features of the oligonucleotide composition are responsible for clade clustering. Additionally, BLSOM also provided information concerning the special genomic region possibly undergoing RNA modifications. CONCLUSIONS: The BLSOM has powerful image display capabilities and enables efficient knowledge discovery about viral evolutionary processes, and it can complement phylogenetic methods based on sequence alignment.

A selective sweep in the Spike gene has driven SARS-CoV-2 human adaptation.

The coronavirus disease 2019 (COVID-19) pandemic underscores the need to better understand animal-to-human transmission of coronaviruses and adaptive evolution within new hosts. We scanned more than 182,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes for selective sweep signatures and found a distinct footprint of positive selection located around a non-synonymous change (A1114G; T372A) within the spike protein receptor-binding domain (RBD), predicted to remove glycosylation and increase binding to human ACE2 (hACE2), the cellular receptor. This change is present in all human SARS-CoV-2 sequences but not in closely related viruses from bats and pangolins. As predicted, T372A RBD bound hACE2 with higher affinity in experimental binding assays. We engineered the reversion mutant (A372T) and found that A372 (wild-type [WT]-SARS-CoV-2) enhanced replication in human lung cells relative to its putative ancestral variant (T372), an effect that was 20 times greater than the well-known D614G mutation. Our findings suggest that this mutation likely contributed to SARS-CoV-2 emergence from animal reservoirs or enabled sustained human-to-human transmission.

Human cell-dependent, directional, time-dependent changes in the mono- and oligonucleotide compositions of SARS-CoV-2 genomes.

BACKGROUND: When a virus that has grown in a nonhuman host starts an epidemic in the human population, human cells may not provide growth conditions ideal for the virus. Therefore, the invasion of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which is usually prevalent in the bat population, into the human population is thought to have necessitated changes in the viral genome for efficient growth in the new environment. In the present study, to understand host-dependent changes in coronavirus genomes, we focused on the mono- and oligonucleotide compositions of SARS-CoV-2 genomes and investigated how these compositions changed time-dependently in the human cellular environment. We also compared the oligonucleotide compositions of SARS-CoV-2 and other coronaviruses prevalent in humans or bats to investigate the causes of changes in the host environment. RESULTS: Time-series analyses of changes in the nucleotide compositions of SARS-CoV-2 genomes revealed a group of mono- and oligonucleotides whose compositions changed in a common direction for all clades, even though viruses belonging to different clades should evolve independently. Interestingly, the compositions of these oligonucleotides changed towards those of coronaviruses that have been prevalent in humans for a long period and away from those of bat coronaviruses. CONCLUSIONS: Clade-independent, time-dependent changes are thought to have biological significance and should relate to viral adaptation to a new host environment, providing important clues for understanding viral host adaptation mechanisms.

Variable routes to genomic and host adaptation among coronaviruses.

Natural selection operating on the genomes of viral pathogens in different host species strongly contributes to adaptation facilitating host colonization. Here, we analyse, quantify and compare viral adaptation in genomic sequence data derived from seven zoonotic events in the Coronaviridae family among primary, intermediate and human hosts. Rates of nonsynonymous (dN ) and synonymous (dS ) changes on specific amino acid positions were quantified for each open reading frame (ORF). Purifying selection accounted for 77% of all sites under selection. Diversifying selection was most frequently observed in viruses infecting the primary hosts of each virus and predominantly occurred in the orf1ab genomic region. Within all four intermediate hosts, diversifying selection on the spike gene was observed either solitarily or in combination with orf1ab and other genes. Consistent with previous evidence, pervasive diversifying selection on coronavirus spike genes corroborates the role this protein plays in host cellular entry, adaptation to new hosts and evasion of host cellular immune responses. Structural modelling of spike proteins identified a significantly higher proportion of sites for SARS-CoV-2 under positive selection in close proximity to sites of glycosylation relative to the other coronaviruses. Among human coronaviruses, there was a significant inverse correlation between the number of sites under positive selection and the estimated years since the virus was introduced into the human population. Abundant diversifying selection observed in SARS-CoV-2 suggests the virus remains in the adaptive phase of the host switch, typical of recent host switches. A mechanistic understanding of where, when and how genomic adaptation occurs in coronaviruses following a host shift is crucial for vaccine design, public health responses and predicting future pandemics.

Emerging of a SARS-CoV-2 viral strain with a deletion in nsp1.

BACKGROUND: The new Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), which was first detected in Wuhan (China) in December of 2019 is responsible for the current global pandemic. Phylogenetic analysis revealed that it is similar to other betacoronaviruses, such as SARS-CoV and Middle-Eastern Respiratory Syndrome, MERS-CoV. Its genome is ∼ 30 kb in length and contains two large overlapping polyproteins, ORF1a and ORF1ab that encode for several structural and non-structural proteins. The non-structural protein 1 (nsp1) is arguably the most important pathogenic determinant, and previous studies on SARS-CoV indicate that it is both involved in viral replication and hampering the innate immune system response. Detailed experiments of site-specific mutagenesis and in vitro reconstitution studies determined that the mechanisms of action are mediated by (a) the presence of specific amino acid residues of nsp1 and (b) the interaction between the protein and the host's small ribosomal unit. In fact, substitution of certain amino acids resulted in reduction of its negative effects. METHODS: A total of 17,928 genome sequences were obtained from the GISAID database (December 2019 to July 2020) from patients infected by SARS-CoV-2 from different areas around the world. Genomes alignment was performed using MAFFT (REFF) and the nsp1 genomic regions were identified using BioEdit and verified using BLAST. Nsp1 protein of SARS-CoV-2 with and without deletion have been subsequently modelled using I-TASSER. RESULTS: We identified SARS-CoV-2 genome sequences, from several Countries, carrying a previously unknown deletion of 9 nucleotides in position 686-694, corresponding to the AA position 241-243 (KSF). This deletion was found in different geographical areas. Structural prediction modelling suggests an effect on the C-terminal tail structure. CONCLUSIONS: Modelling analysis of a newly identified deletion of 3 amino acids (KSF) of SARS-CoV-2 nsp1 suggests that this deletion could affect the structure of the C-terminal region of the protein, important for regulation of viral replication and negative effect on host's gene expression. In addition, substitution of the two amino acids (KS) from nsp1 of SARS-CoV was previously reported to revert loss of interferon-alpha expression. The deletion that we describe indicates that SARS-CoV-2 is undergoing profound genomic changes. It is important to: (i) confirm the spreading of this particular viral strain, and potentially of strains with other deletions in the nsp1 protein, both in the population of asymptomatic and pauci-symptomatic subjects, and (ii) correlate these changes in nsp1 with potential decreased viral pathogenicity.